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ABSTRACT 

A variation of Zeilberger's holonomic ansatz for symbolic de- 
terminant evaluations is proposed which is tailored to deal 
with Pfaffians. The method is also applicable to determi- 
nants of skew-symmetric matrices, for which the original 
approach does not work. As Zeilberger's approach is based 
on the Laplace expansion (cofactor expansion) of the deter- 
minant, we derive our approach from the cofactor expansion 
of the Pfaffian. To demonstrate the power of our method, we 
prove, using computer algebra algorithms, some conjectures 
proposed in the paper "Pfaffian decomposition and a Pfaffian 
analogue of q-Catalan Hankel determinants" by Ishikawa, 
Tagawa, and Zeng. A minor summation formula related to 
partitions and Motzkin paths follows as a corollary. 

Categories and Subject Descriptors 

G.2.1 [Discrete Mathematics]: Combinatorics — Recur- 
rences and difference equations; G.4 [Mathematical Soft- 
ware]: Algorithm design and analysis 

General Terms 

Algorithms, Theory 

Keywords 

Pfaffian, determinant, minor, holonomic systems approach, 
WZ theory, symbolic summation, computer proof, Motzkin 
number 

1. INTRODUCTION 

Pfaffians are a very important concept in combinatorics 
and in physics, for example, for the enumeration of plane 
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partitions, Kasteleyn's method for the dimer models, etc. 
We introduce an algorithmic method for evaluating Pfaffi- 
ans which allows us to solve such problems automatically by 
computer; we demonstrate its applicability by proving a few 
conjectures in [6], concerning Pfaffians of interesting combi- 
natorial numbers. Our approach is a variation of Zeilberger's 
holonomic ansatz for evaluating determinants, which we re- 
call in the following, for sake of self-containedness. 

In |17| , Zeilberger proposed an algorithmic approach for 
evaluating and/or producing rigorous proofs of determinant 
evaluations of the form 

det (fflij)i<i,j<n = b n 

(we use checked letters here to avoid confusion with the 
quantities introduced in Section [2]). The goal is achieved in 
a completely automatic fashion, using computer algebra al- 
gorithms for guessing recurrences and symbolic summation. 
The key point is to guess [8] a suitable (implicit) description 
of an auxiliary function c„,j and then prove that it satisfies 
the three identities 

&n,n = 1 (n > 1), (1) 

n 

} J Cn,ja,i,j =0 (1 < i < n), (2) 

i=i 

" b 
^Cn,ja n , 3 = j-^— (n ^ 1). (3) 

The determinant evaluation follows as a consequence, using 
Laplace expansion w.r.t. the last row and induction on n. 

In principle, the approach is applicable if the matrix is 
never singular, i.e., if b n £ for all n ^ 0. But in or- 
der to turn Identities H}-© m t° routinely provable tasks, 
Zeilberger additionally requires that the matrix entries ht,j 
constitute a bivariate holonomic sequence and that the ra- 
tios of two consecutive determinants b n /b n -i form a uni- 
variate holonomic (P-finite) sequence (in other words, b n is 
required to be what one could call hyper-holonomic) . This 
is the reason why he termed his approach the holonomic 
ansatz |17l I16| . But still, even if all these conditions are 
satisfied, Zeilberger's holonomic ansatz is not guaranteed to 
succeed, because it relies on the fact that the auxiliary func- 
tion Cnj (that is, the cofactors of the Laplace expansion 
with respect to the last row of the n x n matrix, divided by 
the determinant of the (n — 1) X (n — 1) matrix) turns out 
to be holonomic, too. This may be the case or not. If one 
is lucky, i.e., if c„j satisfies sufficiently many linear recur- 
rence equations with polynomial coefficients and therefore 



is holonomic, then the holonomic machinery will produce a 
P-finite recurrence for the sum on the left-hand side of Q. 
Such a recurrence can then be used to prove a (conjectured) 
determinant evaluation b n by substituting the ratio b n /b n -i 
into this recurrence and comparing initial values, or even, if 
the recurrence is not too complicated, to solve it explicitly 
and obtain a closed form for the determinant. 

For a more detailed description and justification of the 
holonomic ansatz, see |17l 112] . We also recommend the 
beautiful essay [14] for the reader who is interested in deter- 
minant and Pfaffian evaluations in general. 

In the present paper, we introduce a variation of Zeil- 
berger's method that is tailored particularly for Pfafnans. 
Recall that Pfafnans are defined only for skew-symmetric 
matrices and that the square of the Pfaffian equals the de- 
terminant. As a trivial consequence, our approach addresses 
determinants of skew-symmetric matrices as well. Clearly 
Zeilberger's holonomic ansatz cannot be applied to skew- 
symmetric matrices, since the determinant in this case van- 
ishes whenever the dimension is odd. Another extension 
of Zeilberger's ansatz, the so-called double-step method, is 
applicable to matrices that are zero either for even or odd 
dimensions 13 . Concerning the evaluation of determinants 
only (not Pfafnans), the double-step method is more general, 
as it does not assume skew-symmetry, but at the same time 
much more complicated and less efficient than our approach 
for Pfafnans. 

In the next section, we state the cofactor expansion of 
the Pfaffian and use it to develop our algorithmic approach 
for dealing with evaluations of Pfaffians; this means proof 
and/or discovery, as in Zeilberger's approach for determi- 
nants. In the following Sections [3H5] this method is used to 
solve some open problems posed in [6]. The details of our 
computer proofs are provided as supplementary electronic 
material on the webpage 

http : //www. rise . jku. at /people/ ckoutsch/pf af f ians/ 

in form of a Mathematica notebook. It is supposed to en- 
able the reader to reproduce our results and do further ex- 
periments. In Section [S]we use our results (Theorem [2]) to 
prove an interesting minor summation formula where the 
sum ranges over certain partitions and the matrix entries 
are variations of Motzkin numbers. We conclude this article 
by posing some open problems as future challenges. 

2. PFAFFIANS 

Let n be a positive integer and let A = (ai,j)i^i,j^2n be 
a 2n by 2n skew-symmetric matrix, i.e., dj,i = —cnj, whose 
entries Oj,j are in a commutative ring. Note that it is com- 
pletely determined by its upper triangular entries dij for 
1 < i < j < 2n. The Pfaffian Pf (A) of A is defined by' 

Pf(-A) = ^^e((7i,<72, • ■ • , CT2n-l) C2n) «ctict 2 • • • a CT 2 „_ 1 <T 2 , l - 

where the summation is over all partitions 
{{ffl, er 2 }, . -. , {cT2n-l,CT2n}} 

of [2n] — {1,2,..., 2n} into two-elements subsets, and where 
e(ci, (72, ■ ■ ■ , C2n-i) o~2n) denotes the sign of the permutation 

( 1 2 • ■ • 2n - 1 2n \ 

\ai (72 ■■■ (?2n-l <T2n J 



A permutation (<ti,<T2, • ■ • , f2n-i , o" 2n ) which arises from a 
partition of [2n] into 2-elements blocks is called a perfect 
matching or a 1-factor. For any permutation n of [2n], let 
A w = {a^uwtj)) denote the skew-symmetric matrix obtained 
by the natural action of n on both rows and columns. From 
the definition above it is easy to see that 

Pf(yT) = sgnvr Pf(A). 

Hence, if any two rows and/or columns are coinciding in A, 
the Pfaffian Pf(^4) of A vanishes. It is a well-known fact that 
Pf(^4) 2 = det(A). Now let I = {h, . . . , i r } be an r-element 
subset of [2n]: we denote by 

A(I) = A{i 1 ,...,i r ) 

the skew-symmetric (2n — r) x (2n — r) matrix obtained 
from A by removing the rows ii,...,i r and the columns 
ii, . . . , i r . Also let us define Tij for 1 ^ i, j ^ 2n by 

f(_l)J-*-ipf A (i,j) ifi<j, 
r w = ^(-l) i -»'PfA(j,i) ifj<i, 
(o ifi = j. 

The Laplace expansion formula for Pfaffians reads as follows. 

Proposition 1. Let A — (fli^)i^i,i^tn be a skew-sym- 
metric matrix, and Tij be as above. Then we have 

2n 2n 

/] a it k^j,k = 2^ ak,iTk,j — Sij Pf A. 
fe=i fe=i 

Proof. This statement and its proof are found in [7]. □ 

Hence if one puts b^n — Pf A = Pf (aj,j)i^i,j^2n and 
/F 2n _i,2n for 1 < i < 2n— 1, then Proposition [T] 
implies that C2n,j satisfies the following three identities: 

C2n,2n-1 =1 (ll > 1), (1) 

2n— 1 

cin^ai.j =0 (1 ^ j < 2n), (2) 

i=l 
2n-l 

c 2 „.iai,2n = 7——— (n > 1). (3) 

r- i 02n-2 
%— 1 

Conversely, one easily sees that the bivariate sequence C2n,j 
is uniquely characterized by Equations ([l]) and ([2]), and we 
can regard H])-® as the formulation analogous to H}-(|!} 
in order to evaluate the Pfaffian Pf A = Pf(oi,j)i<i.j'<2w 
For this purpose, one first has to guess a suitable implicit 
(i.e., holonomic) description of the function C2n,i and then 
show that it indeed satisfies the above identities. Induction 
on n concludes the proof. The methodology is illustrated in 
detail by an example in Section [3] 

Identities fl]), 121), and ((3| can be proven algorithmically 
in the spirit of the holonomic systems approach [16]. In 
the following sections the software package HolonomicFunc- 
tions 11 which runs under the computer algebra system 
Mathematica is employed for carrying out the necessary 
computations. The thesis [9] describes the theoretical back- 
ground and the algorithms implemented therein. 

3. A MOTZKIN NUMBER PFAFFIAN 

This section gives a detailed computer proof of a Pfaf- 
fian involving the Motzkin numbers. It is stated as an open 



problem in |6j, see Formula (6.3) there. The Motzkin num- 
bers M n can be obtained by the formula 

-s^(;)(:)-(-*.^). 

where 2F1 stands for the Gaufi hypergeometric function. 
They count Motzkin paths from (0, 0) to (n, 0); recall that a 
Motzkin path is a path in the lattice N 2 , that uses only the 
steps 

U = (l,l), #=(1,0), D = (l,-1) 
and never runs below the horizontal axis (see [4]). 

Theorem 2. For all integers n 1 the following identity 
holds: 

n-1 

pf (o- - m+i-s) 1<itj<2n = n ( 4fc + !)• ( 4 ) 

fc=0 

Proof. The proof is split into several parts which are 
presented in the Sections I3.1H3.5I below. The details of the 
computations are contained in the supplementary electronic 
material mentioned in the introduction. □ 

3.1 Implicit Description for c 2n ,i 

The first step is to determine the auxiliary function C2n,i 
that appears in identities (jXJ) — <(3j ; where now 

a i,j = U ~ i)M l+j _ 3 . 

Using the method of guessing, as implemented in the Mathe- 
matica package Guess [8], one comes up with an implicit 
description of this unknown function, namely the follow- 
ing three linear recurrence equations with polynomial co- 
efficients: 

(i-l)(2n-3)(4n-7)c 2n ,i = 

-{2n + i - 4)(8in - 8i - 8n 2 + 6n + 3)c 2(n _i),,_i + 
(i - l)(16m - 16i + 8n 2 - 34n + 27)c 2(n _ 1)>i + (5) 
24i(z - l)(n - l)c 2 („_i) li+ i - 
(2n - 3)(4n - 7)(2n - i)c 2n ,i-i, 

(n-2)(2n-5)(4n-ll)(4n-7)(2n— i-2)(2n-i-l)c2n,i = 
(2n-5)(4n-ll)(8i 2 n 2 - 24i 2 n + 17i 2 - 16in 2 + 48m - 

33i - 16ti 4 + 108n, 3 - 258?i 2 + 258™ - 92)c 2(n _ 1) , 4 - 
(n - l)(4n - 7)(2n + i - 5)(32in 2 - 122m + (6) 

117i - 32n 3 + 168n 2 - 280n + 144)c 2 ( n _ 2 ),i + 
6i(4i + l)(n - 2)(n - l)(2n - 3)(4n - 7)c 2(n _ 2) , l+1 + 
36i(i + l)(n - 2)(n - l)(2n - 3)(4n - 7)c 2(n _ 2)jl+2 , 

18n(i -3)(i- 2)(z - l)c 2n ,i = 

(2n + i - 4)(10i 2 n - 24m 2 - 63m + i + 16n 3 + 

76n 2 + 97n - 3)c 2w ,j-3 - 
2(i - 3)n(7i 2 - 12m - 46i + 33n + 73)c 2n ,;_ 2 - 1 ' 

3(i - 3)(i - 2)n(14i - 12n - 39)c 2rM -i - 
(2n - l)(4n - 3)(2n - i + 4)(2n - i + 3)c 2( „ +1)ii _ 3 - 

When they are rewritten in operator notation, these re- 
currences form a left Grobner basis in the corresponding 
noncommutative operator algebra, which is a bivariate poly- 
nomial ring (in the indeterminates Si and S n , denoting the 
forward shift operators w.r.t. i and n, respectively) with co- 
efficients in Q(i, n). Together with the initial values 

C2,l = 1, C 2j2 = C 2j 3 = 0, C4,l = 2, 



they uniquely define the bivariate sequence (o2 n ,i)n,i^i- Note 
that the leading coefficients of ([5]), ([6]), and |[7| never van- 
ish simultaneously in the region where these recurrences are 
used to produce the values c 2n ,i (in the first quadrant, basi- 
cally). 

3.2 Boundary Conditions 

The sequence c is now extended to (c 2n ,i) n >i,iez and it is 
proven that the assumption c 2n ,i = for i ^ and for i ^ 2n 
is compatible with the recurrences (H}-©. This knowledge 
will be useful for the subsequent reasoning. 

Provided with the appropriate initial conditions (c 2 ,o = 
C4.0 = 0), it is obvious that the recurrence produces zeros 
on the line i = 0, since the terms c 2 ( n _ 2 ) ii+ i and c 2 ( n _ 2 ) j i +2 
vanish. Similarly for i = — 1, since the term c 2 ( n _ 2 ) i+2 still 
vanishes; again assuming c 2 ,_i = C4.-1 = 0. Because of 
these two zero rows, it is clear that everything beyond them 
(i.e., for i < — 1) must be zero as well. A simple computation 
shows that setting the initial conditions to is compatible 
with the recurrences ©-([7]). 

Since the leading coefficient of ([5j does not vanish for any 
integer point in the area n 2 and i ^ 2n, this recurrence 
can be used to produce the values of c 2n ,i in this area. The 
support of ([5]) indicates that only c 2 „ i2n = needs to be 
shown. The first instances of this sequence are zero by con- 
struction, and thus we have just to check that the third-order 
recurrence (not printed here) that is automatically derived 
for C2n,2n does not have a singularity in its leading coeffi- 
cient; this is indeed not the case. The univariate sequence 
c 2n ,2n is called the diagonal of the bivariate sequence c 2n ,i. 
Diagonals appear frequently in combinatorial problems and 
their fast computation is a topic of ongoing research in com- 
puter algebra. We used the command DFiniteSubstitute 
of pT| here to perform the substitution i — > 2n, which cor- 
responds to the computation of the diagonal. 

It remains to show that c 2 ,j = for i > 2, which is done 
in a similar fashion. 

3.3 Identity CQ) 

Analogously to the computation of the diagonal in the 
previous section, an annihilating operator for c 2nj2n _i (of 
order 4, not printed here) is obtained. Its leading coefficient 
has no nonnegative integer roots, and it has the operator 
S n — 1 as a right factor. Therefore it annihilates any constant 
sequence. The four initial values are 1 by construction and 
therefore c 2n ,2n-i = 1 for all n 6 N. 

3.4 Identity © 

Once the implicit descriptions of the bivariate sequences 
aij and C2n,i are available, in terms of zero-dimensional left 
ideals of recurrence operators, the summation identities ((JJ 
and ([3]) are routinely provable, thanks to software packages 
like HolonomicFunctions [11] . The strategy is as follows: 
first the closure properties of holonomic functions are em- 
ployed to compute recurrences for the product C2 n ,iaij; the 
command DFiniteTimes does the job. Then the method 
of creative telescoping is invoked to produce some recur- 
rences for the left-hand side of @ (this expression is de- 
noted by g„ t j in the following). Two different algorithms for 
this task are implemented in our package, namely the com- 
mands CreativeTelescoping (Chyzak's algorithm [3]) and 
FindCreativeTelescoping (an alternative ansatz proposed 
by the second author [10]). In order to prove Identity ((2| 



for instance, some operators of the form 

P(j, n, S h S n ) + (Si ~ l)Q(i,j, n, St, S s ,S n ) 

which annihilate the summand C2n,i<Jij are computed. It 
has already been proven in Section[3]2]that C2n,i is zero out- 
side the summation range which implies that the sum runs 
over natural boundaries. Therefore the principal parts (or 
telescopers, denoted by P above) of the creative telescoping 
operators annihilate the sum, and the delta parts (denoted 
by Q) can be disregarded. As a result we find 

j(An-7)(2n + j-2)g n , j = 

j(An - 3)(j - n + l)g n -i,j + 

(n - l)(4n - 3)(2n -j - 3)g n -i, j+ i, 

(j-2n)(2n+j-2)g nJ = 

3(j - 2)(j - l) fl „j_ a + (j - l)(2j - 3)s„j_i. 

A close inspection reveals that only the initial values <?i,i, 
g2,i, and t/2,2 need to be given, if the above recurrences shall 
be used to compute all values of g n j for n ^ 1 and 1 ^ 
j < 2n. A simple calculation shows that they are all zero, 
concluding the proof of 

Note that the above reasoning is somehow about the max- 
imal possible area: if one tries to extend it further, i.e., to 
show that g n j — in the whole first quadrant, the first step 
being the points j — 2n, then the second recurrence, the 
only one that is applicable in this case, breaks down. In- 
deed, the values g n ,2n are nonzero as is demonstrated in the 
next section. 

3.5 Identity © 

Identity ((3} is done in a very similar fashion, using the 
method of creative telescoping. Again the summation is over 
natural boundaries. Thus the principal part of the computed 
creative telescoping operator gives rise to a recurrence for 
the left-hand side of <(3j which is denoted by r n here: 

2(4n - 11) (4n - 7)(4n - 5)(7n - 13)r„ = 

(An - ll)(350n 3 - 1413n 2 + 1798™ - 714)r n _i - 
9(n - 2)(2n - 3)(4n - 7)(7ti - 6)r n _ 2 . 

For n = 1 and n — 2 the summation in ((3| yields the initial 
values n = 1 and r : 2 = 5. It is easily verified that the unique 
solution of the above recurrence is r n = An — 3. Since r n = 
fe2n/fe2n-2 gives the quotients of two consecutive Pfaffians 
fen = Pf(ai,i)i<ij<2ni it follows that 

n . n n—1 

6 -=n^=n( 4fc - 3 )=n( 4fc + 1 )- 

k=l k=l k=0 

This concludes the proof of Theorem 

4. A DELANNOY NUMBER PFAFFIAN 

We now consider a Pfafnan that appears as Formula (6.4) 
in [6], again as an open problem. 

Theorem 3. Let 



denote the n-th central Delannoy number. Then for all in- 
tegers n 1 the following identity holds: 

n—l 

Pf((i-i)A + ,-- 3 ) 1 ^ > , <2n =2(" +1 H«- 1 )(2n-l)n(4fc-l). 

fc=i 

(8) 

Proof. The auxiliary function C2n,i in this example can 
be defined by the following recurrences (plus a sufficient 
amount of initial values): 

2(i-3)(i-2)(i-l)c2n,i = 

3(i - 3)(i - 2)(8i - 27)c 2n ,i_i - 

(i - 3)(76i 2 - 589i - 8n 2 + 16n + 1109)ca»,i-2 + 

3(8i 3 - 105i 2 - 16m 2 + 32m + 443i + 68n 2 - 

136n — 600)c2n,<-3 — 
(2i - 11) (i -2n- 3)(i + 2n~ 7)c 2n ,i-i, 

2(n - 2)(2n - 3)(4n - 9)(i - 2n + l)(i - 2n + 2)c 2 „,i = 
(n - l)(i + 2n - 5)(68i 2 n - 102i 2 - 96m 2 + 178m - 

43i + 64n 3 - 208n 2 + 200n - 56)c 2(n _ 1)ii - 
6i(n-l)(2n - 3)(35i 2 + Ain - 66i - n + 14)c 2 ( re _i), i+1 + 
i(i + l)(n - l)(2n - 3)(70z + An - 31)c 2 ( n _i), l+ 2 -' 
6i(» + l)(t + 2)(n - l)(2n - 3)c 2( „_ 1)]i+3 . 

The proof is very analogous to the one of Theorem[2] see the 
accompanying Mathematica notebook for the details. □ 

5. A NARAYANA NUMBER PFAFFIAN 

The following Pfafnan appears as Formula (6.6) in [6]: 

Theorem 4. Let N„(x) denote the n-th Narayana poly- 
nomial defined by 

N (x) = 1, 

«.« - U^-h <nsl) - 

Then for all n ^ the following identity holds: 

n-l 

pf (cj - m + ^))^ 2n = ^ IK 4fc + !)• (9) 

fc=0 

Proof. Again, the proof of this evaluation is analogous 
to the previous ones of ((4]) and see the accompanying 
Mathematica notebook for the details. The main difference 
is that now the free parameter x is involved, which on the one 
hand makes the computations and the intermediate results 
more voluminous. On the other hand, some arguments in 
the proof (like "the leading coefficient of some recurrence is 
never zero") become more intricate. 

One solution to address the latter issue is to argue that 
a: is a formal parameter; then any polynomial in x which 
is not identically zero, is considered to be nonzero (as an 
element in the corresponding polynomial ring). If one feels 
uneasy about this argument, one can as well try to find 
conditions under which all steps of the proof are sound; for 
our reasoning the assumption x < — 1 was sufficient. Hence 
the evaluation is proven only for x < — 1. But for specific n, 
the Pfaffian is a polynomial in x (of a certain degree), as 
well as the evaluation on the right-hand side of (0). Thus 
their difference is a polynomial in x which has been proven 
to be zero for all a; < —1. By the fundamental theorem of 
algebra it follows that this polynomial is identically zero, 



and therefore the evaluation of the Pfafnan is true for all 
complex numbers x. □ 



Corollary 5. Let 



^-g*ir( B i*) 



n + k \ [2k 
k 



denote the (large) Schroder numbers. Then for all integers 
n ^ the following identity holds: 



2" 2 TJ(4*+1). 



Proof. This identity follows from Theorem [4] and the 
equality S„ = AT n (2); the latter fact can be easily proven 
from the definitions of these quantities using Zeilberger's 
algorithm, for example. □ 

In [6] it has already been noted that Theorem [4] implies 
the Pfafnan of Corollary [S] involving the Schroder numbers. 
Similarly, it is stated there that also the Pfafnan Q is a 
special case of Theorem [4] However, in order to reflect the 
historic evolution of our results and for reasons of a clear 
presentation, we included Theorem [2] and its detailed proof 
in this article. 

6. APPLICATION OF THEOREM 2 

Let A = (a«,j)i^i<w, j>i be any n-rowed matrix. If J = 
{ji, . ..,j n } is a set of column indices, then we write Aj = 
Aj 1: ...j n for the square submatrix of size n obtained from A 
by choosing the columns indexed by J. If A — (aj, 3 )ij>i is 
a matrix with infinitely many rows and columns, and 7 = 
{ii,. . . ,i n } (resp. J = {ji, . . . , j n }) is a set of row (resp. 
column) indices, then let Aj = Aj 1 '"'^ denote the square 
submatrix of size n obtained from A by choosing the rows 7 
and the columns J. 

A partition is a nonincreasing sequence A = (Ai, A2, . . . ) 
of nonnegative integers with only finitely many nonzero el- 
ements. The number of nonzero elements in A is called the 
length of A and is denoted by l{\). An odd partition is a 
partition with odd parts and an even partition is a par- 
tition with even parts. The conjugate of A is the parti- 
tion A' = (A'x, A2, . . . ), where A; is the number defined by 
\' t — I Aj ^ i}. Given a partition A such that l(X) ^ n, 
let 7„(A) denote the n-element set of nonnegative integers 
defined by 

7„(A) = {A n + 1, A„_i + 2, . . . , Ai + n} . 

For example, A = (3, 3, 1, 1) is an odd partition of length 4, 
and 7 4 (A) = {2, 3, 6, 7}. The conjugate of A equals (4, 2, 2), 
which is an even partition. 

Let 77(n) = (h(i,j)) 1<i<n j >1 denote the n-rowed matrix 
whose entries are given by 



fc (i,»-l)=( fc _ 1 ) a f 1 ( — ;4 

/ •_ o\ / k-i+l k-i+2 

Ki,2k) = (i-1)(-_ 1 ) 2 F 1 ( ;4 



(10) 
(11) 



In fact h(i, 2k — 1) is the number of Motzkin paths from (0, 0) 
to (i - 1, k - 1). We also note that h(i, 2k) = fc[z i+fe - 1 ](l + 
x + x 2 ) 1-1 , where [x n ]f(x) denotes the coefficient of x n in a 



polynomial f(x). For example, if n — 4, then we have 

/10000000..A 
, , _ 1 1 1 0.. 
H( - i >~ 22221000.. 

\4 6 5 6 3 3 1 .../ 

For example, h(A, 3) = 5 gives the number of Motzkin paths 
from (0,0) to (3,1): 

UUD, (77777, UDU, HUH, HHU. 

Meanwhile, h(4, 4) = 6 equals 2 times the coefficient of x 5 in 
(1 + x + x 2 ) 3 . The main purpose of this section is to give a 
proof of the following theorem as a corollary of Theorem [5] 

Theorem 6. Let n be a positive integer, and let H(n) be 
as above. Then we have 



det77(2n) /2 „ (A) 

A 

A, A 7 even 



n( 4fc+i )> 



(12) 



where the sum on the left-hand side runs over all even par- 
titions A such that £(X) ^ 2n and A' is also even. 

We notice that this theorem is a consequence of an addition 
formula for 27^1 and Theorem [2] But it is not so easy to 
find a lattice path interpretation of det H(2n)i 2n r\-\ since 
we do not know a lattice path interpretation of h(i, 2k). To 
prove the theorem, we cite the following two lemmas from [5] 
and 0. 

Lemma 7. If i and j are nonnegative integers, then we 
have 



E 



1F1 



k—i+l k—i 
2 ' 2 .4 
k + 2 ' 



1F1 



k—j+1 k—j 
2 ' 2 

k + 2 



1 — i— j 
2 ' 



;4 



(13) 



Proof. The proof can be found in [S] Lemma 5.2]. □ 

Lemma 8. For n £ N let T — (fij)i^j^2n, j^i be an n- 
rowed matrix, and let A — (a<,j)j,j>i be a skew-symmetric 
matrix with infinitely many rows and columns, i. e. aj,i ~ 
—di t j for i, j > 1. Then we have 



^2 Pf(Af)det(I>) = Pf(0), 



(14) 



where the sum on the left-hand side runs over all 2n-element 
sets of positive integers and the skew- symmetric matrix Q 
is defined by Q = (Qij) = TAT T whose entries may be 
written in the form 



Qi,j = 52 a k ,i det(T£ 



), (l^ij^n). 



(15) 



Proof. The proof of this minor summation formula can 
be found in [7] Theorem 3.2]. □ 

Proof of Theorem El Set 

1 if i = 2k — 1 and j — 2k for some k G N, 
a», j = \ — 1 if i = 2k and j — 2k — 1 for some t£N, 
otherwise. 



and Uj = h(i,j) in (fl4)) where h(i,j) is as denned in (flO]) 
and (|11[) . Then one can show by direct calculation Qi t] - 
in (|15[) is given by 



^-2- det U(j,2fe-l) h(j,2fc) 



o->E(::i)fi:r 



fc-j+l fc-j+2 

fc + 1 

/ — — L!- - — — 

/.• i y • ' \ /.■ - i 



x a *| 2 fc+ 2 ;4 



By (|13p we obtain 



Qij = <j-fyFA 2 ' 2 ;4) =(j-i)M i+j _ 3 



3 — i — j 4— <— jj 
2 ' 2 
2 



which, using (|14[1 . gives 

^ Pf (^) det#(2n), = Pf((j - OM^-s) 



It is not hard to see that Pf(Aj) = 1 if I = /2u(A) for a 
partition A such that £(A) ^ 2n and A, A' are even, and 
Pf(-Aj) = otherwise (see [7]). Hence we obtain the desired 
formula (|12p as a consequence of Theorem [2] □ 



We can regard the numbers 



h(i,2k - 1) = 



z - 1 
k- 1 



2F1 



fc — i fc — B-j-l 



fc + 1 



2 .4 



as a generalization of the Motzkin numbers M n since they 
count the Motzkin paths from (0,0) to (i — 1, k — 1), and 
write = h(i, 2k — 1) hereafter. In fact (A^ ) gives 

the (k — l)-th column of the Motzkin triangle [4]. Note that 
M n = so that Theorem [2] reads 

n-l 

pf((i-i)x^ f _ 2 ) =n(4*+i). 

It may now be attractive to present a generalization of The- 
orem [5] as follows. 

Conjecture 9. Let n and k be positive integers. 

(i) Then the Pfaffian 

Pf(jj-i)M 

equals 

m — 1 k — 1 

n Y[{Aki+2j+k) 

i=0 j=0 

if m = n/k is an integer, and it equals 

j=l J / \ i=0 3=1 

if k is odd and m — (n+ L&/2J) Ik is an integer. The 
Pfaffian is zero in all other cases. 



,(fe) 
H+3-2 



l^i,3^2n 



(ii) Meanwhile, the Pfaffian 
equals 

m — 1 k — 1 

n y[(ah + 2j + k + 1) 

i=0 j=0 

if m = rz/fc is an integer, and it equals 

k/2 \ / ra-1 k 



i=0 j = l 



z/fc is euen ana! m = (n + k/2) Ik is an integer. The 
Pfaffian is zero in all other cases. 

We want to remark that for k — 1 part (i) is just Theorem [5] 
and part (ii) can be proven analogously. Unfortunately these 
two Pfafnans are periodically zero if k ^ 2. This prevents us 
from applying our method to the conjecture, since we con- 
sider the quotient of two consecutive Pfafnans. Of course, 
one could come up with a Pfaffian analogue of the double- 
step method presented in |13| . which would settle Conjec- 
ture [§] for the special case k = 2. This construction may 
be extended for k — 3, k = 4, etc., at the cost of more and 
more involved computations. However, this approach will 
not work for symbolic k in general. 

Conjecture [9] can be regarded as a Pfaffian analogue of 
the Hankel determinants of Motzkin numers [T] Proposi- 
tion 2] , and Hankel determinants of sums of two consecutive 
Motzkin numbers [5J Theorem 3.2]. Many combinatorial ar- 
guments are known for the Hankel determinants, but little 
is known for Hankel Pfafnans (see [HI US]). It may be inter- 
esting to discover a combinatorial reason why we can expect 
such a nice formula for the Hankel Pfafnans of (sums of) 
Motzkin numbers. 
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